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Chloroquine,  long  the  default  first-line  treatment  against  malaria,  is  now  abandoned  in  large  parts  of  the 
world  because  of  widespread  drug-resistance  in  Plasmodium  falciparum.  In  spite  of  its  importance  as  a 
cost-effective  and  efficient  drug,  a  coherent  understanding  of  the  cellular  mechanisms  affected  by 
chloroquine  and  how  they  influence  the  fitness  and  survival  of  the  parasite  remains  elusive.  Here,  we 
used  a  systems  biology  approach  to  integrate  genome-scale  transcriptomics  to  map  out  the  effects  of 
chloroquine,  identify  targeted  metabolic  pathways,  and  translate  these  findings  into  mechanistic  in¬ 
sights.  Specifically,  we  first  developed  a  method  that  integrates  transcriptomic  and  metabolomic  data, 
which  we  independently  validated  against  a  recently  published  set  of  such  data  for  Krebs-cycle  mutants 
of  P.  falciparum.  We  then  used  the  method  to  calculate  the  effect  of  chloroquine  treatment  on  the 
metabolic  flux  profiles  of  P.  falciparum  during  the  intraerythrocytic  developmental  cycle.  The  model 
predicted  dose-dependent  inhibition  of  DNA  replication,  in  agreement  with  earlier  experimental  results 
for  both  drug-sensitive  and  drug-resistant  P.  falciparum  strains.  Our  simulations  also  corroborated 
experimental  findings  that  suggest  differences  in  chloroquine  sensitivity  between  ring-  and  schizont- 
stage  P.  falciparum.  Our  analysis  also  suggests  that  metabolic  fluxes  that  govern  reduced  thioredoxin 
and  phosphoenolpyruvate  synthesis  are  significantly  decreased  and  are  pivotal  to  chloroquine-based 
inhibition  of  P.  falciparum  DNA  replication.  The  consequences  of  impaired  phosphoenolpyruvate  syn¬ 
thesis  and  redox  metabolism  are  reduced  carbon  fixation  and  increased  oxidative  stress,  respectively, 
both  of  which  eventually  facilitate  killing  of  the  parasite.  Our  analysis  suggests  that  a  combination  of 
chloroquine  (or  an  analogue)  and  another  drug,  which  inhibits  carbon  fixation  and/or  increases  oxidative 
stress,  should  increase  the  clearance  of  P.  falciparum  from  the  host  system. 

©  2017  The  Authors.  Published  by  Elsevier  Ltd  on  behalf  of  Australian  Society  for  Parasitology.  This  is  an 
open  access  article  under  the  CC  BY-NC-ND  license  (http://creativecommons.Org/licenses/by-nc-nd/4.0/). 


1.  Introduction 

Malaria,  which  is  caused  by  a  protozoan  parasite  of  the  Plas¬ 
modium  genus,  represents  a  serious  global  health  concern  given 
that  nearly  half  of  the  world  population  is  at  risk  of  infection  (WHO, 
2014).  Although  several  anti-malarial  drugs  for  treating  the 
symptomatic  stage  (or  blood  stage)  of  malarial  infection  are 
commercially  available  (Antony  and  Parija,  2016),  their  efficacy  has 
declined  appreciably  in  the  last  few  decades  owing  to  widespread 
drug  resistance  developed  by  the  parasite  (Breman  et  al.,  2004). 
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Plasmodium  falciparum— the  most  lethal  species  of  malar¬ 
ia-causes  approximately  50%  of  all  malarial  infections.  Chloro¬ 
quine  was  the  first-line  malaria  treatment  for  many  decades  until 
drug-resistant  P.  falciparum  strains  became  common.  The  drug 
causes  a  dose-dependent  decrease  in  hemozoin  formation  (Chou 
and  Fitch,  1992;  Slater  and  Cerami,  1992)  and  an  associated  in¬ 
crease  in  toxic  free  heme  in  the  food  vacuole  of  the  parasite 
(Combrinck  et  al.,  2013;  Loria  et  al.,  1999).  Over  the  past  few  de¬ 
cades,  researchers  have  proposed  many  different  mechanisms  for 
chloroquine  action,  including  l)  DNA  intercalation  (Meshnick, 
1990),  2)  alteration  of  digestive  food  vacuole  pH  (Yayon  et  al., 
1985),  3)  inhibition  of  heme  polymerase  (Yayon  et  al.,  1985),  and 
4)  formation  of  a  toxic  chloroquine-ferriprotoporphyrin  IX  complex 
(Sugioka  et  al.,  1987).  Yet,  there  is  no  consensus  on  the  exact 
mechanisms  by  which  chloroquine  kills  the  parasite  and  under 
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what  circumstances  they  operate.  Shedding  light  on  the  possible 
actions  of  chloroquine  is  critical  for  developing  more  potent  drugs 
or  combination  drug  treatments  against  the  parasite  and  for  un¬ 
derstanding  how  the  parasite  can  develop  resistance  against  them. 

Here  we  used  systems  biology  model  and  transcriptional  pro¬ 
files  of  P.  falciparum  to  obtain  information  about  metabolic  and 
biological  precursors  that  determine  the  physiological  or  patho¬ 
physiological  state  of  the  parasite.  Specifically,  we  used  the 
P.  falciparum  transcriptomic  data  obtained  by  Hu  and  colleagues 
during  the  intraerythrocytic  development  cycle  (IDC)  at  different 
concentrations  of  chloroquine  (Hu  et  al.,  2010).  Our  primary  goal 
was  to  use  these  data  to  investigate  the  effect  of  chloroquine  on 
P.  falciparum  DNA  replication.  We  chose  to  study  this  effect  because 
it  has  been  proposed  as  a  mechanism  of  chloroquine  action  in 
different  Plasmodium  species  (Gutteridge  et  al.,  1972;  Meshnick, 
1990;  Polet  and  Barr,  1968a,  b),  in  bacteria  (Ciak  and  Hahn,  1966), 
and  in  bioassays  examining  the  effect  of  chloroquine  on  RNA  and 
DNA  polymerases  (Cohen  and  Yielding,  1965;  O'Brien  et  al.,  1966b; 
Whichard  et  al.,  1972). 

We  first  developed  a  method  that  integrates  the  information 
embedded  in  the  transcriptome  data  with  the  whole-genome 
metabolic  network  model  for  P.  falciparum  (Fang  et  al.,  2014)  to 
predict  the  metabolic  phenotype  corresponding  to  those  tran¬ 
scription  data.  We  validated  the  developed  method  by  using  a  data 
set  from  an  independent  study  that  jointly  capture  transcriptomic 
and  metabolomic  data  from  P.  falciparum,  to  correctly  predict  the 
metabolic  phenotype  tied  to  genetic  perturbation  of  two  Krebs 
cycle  enzymes  (Ke  et  al.,  2015).  We  then  focused  on  the  inhibition  of 
DNA  replication  in  P.  falciparum  (as  a  consequence  of  chloroquine 
treatment)  a  critical  metabolic  phenotype  (Ciak  and  Hahn,  1966; 
Cohen  and  Yielding,  1965;  Polet  and  Barr,  1968a),  and  used  our 
method  to  simulate  the  effect  of  chloroquine  on  the  metabolic 
fluxes  of  P.  falciparum  during  the  IDC.  We  identified  genes  that  were 
substantially  altered  in  response  to  chloroquine  treatment  and 
linked  to  the  inhibition  of  P.  falciparum  DNA  replication.  The  cohort 
of  identified  genes  suggests  that  DNA  replication  is  inhibited  by  the 
downstream  effect  of  heme  accumulation.  Specifically,  our  analysis 
suggests  that  continuous  accumulation  of  heme  inhibits  redox 
metabolism,  carbon  fixation,  and  pyrimidine  metabolism,  which 
leads  to  inhibition  of  DNA  replication  and  facilitates  death  of  the 
parasite.  Our  results  provide  a  mechanistic  explanation  for  why 
parasites  with  an  efficient  redox  metabolism  may  have  a  lower 
sensitivity  to  chloroquine  (Kasozi  et  al.,  2013;  Lehane  et  al.,  2012; 
Meierjohann  et  al.,  2002). 

2.  Materials  and  methods 

2.1.  Transcriptome  data  used  for  simulating  metabolic  flux  profile 
during  IDC 

Metabolic  flux  profiles  of  P.  falciparum  strains  3D7  (Pf3D7)  and 
Dd2  (PfDd2)  during  the  IDC  were  estimated  by  using  transcriptome 
data  previously  reported  by  Llinas  and  colleagues  (Llinas  et  al., 
2006),  who  extracted  cDNA  from  parasites  cultured  in  erythro¬ 
cytes  and  hybridized  them  to  DNA  microarrays  to  yield  hourly 
levels  of  gene  expression  for  the  Pf3D7  (or  PfDd2)  strain  during  the 
IDC.  The  temporal  resolution  of  these  hourly  gene  expression  pro¬ 
files  allowed  us  to  make  high-resolution  time-dependent  metabolic 
flux  predictions  during  the  IDC. 

We  used  previously  reported  gene  expression  data  (Hu  et  al., 
2010)  to  simulate  the  effect  of  chloroquine  on  the  time- 
dependent  metabolic  fluxes  of  Pf3D7  and  PfDd2  during  the  IDC. 
Briefly,  Hu  and  colleagues  (Hu  et  al.,  2010)  administered  chloro¬ 
quine  18  h  after  erythrocyte  invasion  at  concentrations  of  41,  72, 
and  144  nM  for  Pf3D7  and  at  43  nM  for  PfDd2,  and  obtained  gene 


expression  profiles  every  hour  for  8  h  during  the  IDC  following  the 
chloroquine  treatment. 

2.2.  Metabolic  network  model  and  data  processing 

The  metabolic  network  model  used  involved  1045  reactions  and 
376  genes.  The  model  was  identical  to  a  previously  published 
model  (Fang  et  al.,  2014),  with  one  minor  modification:  we  added 
reactions  explicitly  carried  out  by  Pdx2  (PF11_0169)  and  Pdxl 
(MAL6P1.215).  This  was  done  by  1)  modifying  an  old  reaction  that 
incorrectly  assumed  that  pyridoxal  5-phosphate  (PLP)  was  syn¬ 
thesized  from  ribulose  5-phosphate  and  2)  adding  a  new  (and 
separate)  reaction  for  Pdxl  because  experiments  suggest  that  it  can 
synthesize  PLP  independent  of  the  ammonia  generated  by  Pdx2 
(Hanes  et  al.,  2008).  In  the  updated  model,  Pdx2  hydrolyzes  L- 
glutamine  to  yield  glutamate  and  ammonia  while  Pdxl  in¬ 
corporates  the  ammonia  generated  by  Pdx2  to  yield  PLP  from  the 
substrates  ribose  5-phosphate  and  glyceraldehyde  3-phosphate 
(Gengenbacher  et  al.,  2006),  with  the  isomerization  of  ribose  5- 
phosphate  to  ribulose  5-phosphate  as  a  side  reaction. 

We  used  gene-to-reaction  mappings  available  in  the  metabolic 
network  model  to  obtain  reaction  expressions  rexp,  which  depend 
on  the  expression  of  a  gene  (or  genes)  catalyzing  a  particular 
metabolic  reaction.  For  a  reaction  catalyzed  by  a  single  enzyme,  rexp 
was  equal  to  the  expression  of  the  gene  catalyzing  that  reaction. 
However,  if  a  reaction  was  associated  with  more  than  one  gene, 
then  their  gene-to-reaction  mappings  were  defined  by  using  the 
Boolean  operators  “AND”  and  “OR.”  We  implemented  these  oper¬ 
ations  by  taking  the  minimal  (for  the  AND  operator)  and  maximal 
(for  the  OR  operator)  values  of  the  associated  gene  expression  data 
(Song  et  al.,  2014). 

2.3.  Estimating  metabolic  flux  changes  in  response  to  external  stress 

We  recently  developed  a  method  (Fang  et  al.,  2014)  that  in¬ 
tegrates  hourly  gene  expression  data  with  a  whole-genome  meta¬ 
bolic  network  model  to  estimate  the  hourly  metabolic  flux  profiles 
of  P.  falciparum  during  the  48-h  long  IDC  period.  To  take  advantage 
of  additional  studies  that  collect  gene  expression  data  as  a  function 
of  genetic  or  chemical  perturbations,  we  modified  the  method 
because  such  studies  usually  do  not  provide  the  same  granularity 
and  density  of  gene  expression  data  to  fully  cover  the  entire  IDC. 
Therefore,  we  aimed  to  develop  a  method  that  could  utilize  a 
reduced  transcriptomic  data  set  to  predict  alterations  in  metabolic 
flux  profiles  during  the  IDC.  Briefly,  We  used  the  following  steps  to 
estimate  metabolic  fluxes  under  no-stress  conditions: 

Step  1:  Estimate  nutrient  fluxes  sufficient  for  the  nominal 
growth  rate  (40%  of  the  maximum  growth  rate). 

Step  2:  Estimate  reaction  fluxes  for  metabolite  i  corre¬ 
sponding  to  the  nominal  growth  rate  constrained  by  the 
nutrient  fluxes  obtained  in  Step  I. 

Step  3:  Estimate  reaction  fluxes  v\  for  reaction  i  at  any  time  point 
t  during  the  IDC  by  incorporating  time-dependent  tran¬ 
scriptomic  profiles  under  stoichiometric  and  nutrient  uptake 
constraints. 

The  last  step  was  achieved  by  globally  minimizing  for  all  re¬ 
actions  and  time  points  the  value  of  t  that  minimized  the  weighted 

absolute  difference,  XI  \v\  -tfvj  L  where  G  represents  the  set  of 

ieG 

unidirectional  reactions  with  known  gene-to-reaction  mapping. 
The  factor  f\  =  —^4  mm^p;d  3  where,  rlvn ,  represents  reaction 
expression  of  ith  reaction  at  time  t  obtained  from  the  gene-to- 
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reaction  mapping  of  the  P.  falciparum  metabolic  network,  according 
to  the  approach  outlined  in  Section  2.2.  This  factor,  r •,  ensured  that 
the  mean  of  the  ith  reaction  during  the  IDC  was  equal  to  i^.  The 
complete  procedure  has  been  detailed  elsewhere  (Fang  et  al.,  2014). 

In  this  work,  we  estimated  the  organism-level  impact  of  a  drug 
on  metabolism  by  using  transcriptomic  data  obtained  after  drug 
treatment.  The  original  formulation  outlined  above  ( Steps  1—3) 
cannot  be  directly  translated  to  stress  conditions  unless  the  com¬ 
plete  transcriptomic  profile  for  each  stress  condition  during  the 
entire  IDC  is  available  (Fang  et  al.,  2014).  Therefore,  we  introduced  a 
parameter  to  account  for  up-  or  down-regulation  of  individual 
reactions  in  response  to  stress  to  simulate  the  effect  of  exogenous 
stress  on  metabolic  flux  profiles  during  the  IDC.  We  modified  the 

factor  f\  to  where  a,-  denotes  the  median  relative 

change  in  the  reaction  expression  of  the  ith  reaction  following 
chloroquine  treatment,  as  observed  in  previous  experiments  (Hu 
et  al.,  2010).  Note  that  the  mean  of  each  individual  modifying  fac¬ 
tor  r[,  i.e.,  obtained  by  summing  the  values  across  all  time 

points  in  the  IDC,  tracks  with  cq;  that  is,  when  the  mean  is  less  than 
one,  greater  than  one,  or  equal  to  one,  cq<  1,  cq>l,  or  cq  =  1, 
respectively.  In  other  words,  cq  increases  or  decreases  the  nominal 
flux  of  the  ith  reaction  during  the  entire  IDC  depending  on  the 
relative  increase  or  decrease  in  the  ith  reaction  expression.  In 
principle,  a  drug  may  differentially  affect  gene  transcription  at 
different  time  points  during  the  IDC.  Here,  we  assumed  that  the 
median  change  in  gene  transcription  after  a  stress  treatment  was 
representative  of  the  transcriptional  response  throughout  the  IDC. 

2.4.  Validating  whether  reaction  flux  estimates  capture  the  effect  of 
external  stress? 


activated,  because  the  parasite  will  synthesize  them  as  needed 
throughout  the  IDC.  We  used  a  validated  in-house  developed 
method  (Fang  et  al.,  2014)  to  simulate  time-dependent  DNA  syn¬ 
thesis  during  the  IDC.  For  comprehensiveness,  we  introduce  the 
constraint  governing  the  rate  of  DNA  synthesis  in  the  metabolic 
network  model: 

CdATpdATP  +  cdcrPdCTP  +  cdQjpdGTP  +  cdjjpdTTP  —>  DNA  +  (c^^p 

+  cdcrp  +  cdGTP  +  cdTrP)  Diphosphate. 

Here,  dATP ;  dCTP,  dGTP ;  and  dTTP  refer  to  deoxyadenosine, 
deoxycytidine,  deoxyguanosine,  and  deoxythymidine  tri¬ 
phosphates,  respectively.  We  refer  the  reader  to  the  original  article 
(Fang  et  al.,  2014)  for  further  details  on  this  method  for  simulating 
time-dependent  P.  falciparum  macro-molecules  during  the  IDC. 

2.6.  Identifying  genes  critical  for  P.  falciparum  DNA  synthesis 

The  criticality  of  a  gene  is  determined  by  optimizing  the 
following  objective  function:  XI  C/Mp  subjected  to:  SmxnVrva  =0 

i^DNA 

and  vLB  <  v  <  Dub.  Here,  DNA  is  the  gene  set  containing  the  me¬ 
tabolites  that  contribute  to  DNA  synthesis  (introduced  above)  and 
q  is  the  contribution  of  the  ith  metabolite  (M,-).  In  addition,  S  is  the 
stoichiometry  matrix,  v  is  a  vector  containing  all  metabolic  re¬ 
actions  of  the  P.  falciparum  metabolic  network  model,  m  is  the 
number  of  metabolites  in  the  network,  n  is  the  number  of  reactions 
in  the  network,  and  vLB  and  vUB  are  the  upper  and  lower  bounds, 
respectively,  of  metabolic  reactions.  We  deemed  a  reaction  as 
critical  if  there  was  no  DNA  synthesis  after  knocking  out  the  reac¬ 
tion  in  silico.  We  denoted  any  gene  encoding  such  a  reaction  as  a 
critical  gene. 


We  used  a  previously  published  data  set  (Ke  et  al.,  2015)  to  vali¬ 
date  our  methodology  to  simulate  the  effect  of  external  stress  on 
metabolic  fluxes  during  the  IDC.  Ke  and  colleagues  (Ke  et  al.,  2015) 
constructed  a  P.  falciparum  mutant  with  two  tricarboxylic  acid  (TCA) 
cycle  enzymes  deleted  ( AIDH/AKDH ),  cultured  wild-type  and  mutant 
parasites  in  erythrocytes,  and  obtained  matched  transcriptomic  and 
metabolomic  data.  We  used  this  data  set  because  it  describes  changes 
in  gene  expression  (transcriptomics)  and  phenotype  (metabolomics) 
under  the  same  experimental  stress  conditions. 

If  the  reaction  fluxes  estimated  from  previously  reported  tran- 
scriptome  data  (Ke  et  al.,  2015)  adequately  represent  the  reaction 
fluxes  in  vitro ,  then  the  predicted  changes  in  reaction  fluxes  should  be 
sufficient  to  predict  the  changes  observed  in  the  corresponding 
metabolomic  data.  We  used  a  computer  model  of  mitochondrial 
oxidative  phosphorylation  and  the  TCA  cycle  (Wu  et  al.,  2007)  to 
predict  changes  in  TCA  cycle  metabolites  in  response  to  changes  in 
these  estimated  reaction  fluxes  (Ke  et  al.,  2015).  The  kinetic  param¬ 
eters  of  this  computer  model  were  determined  and  validated  by 
using  extensive  experimental  data.  We  simulated  wild-type  re¬ 
actions  by  running  this  model,  with  its  default  parameters,  for  2000  s 
(time  required  to  reach  steady-state  condition)  and  stored  the  con¬ 
centrations  of  TCA  cycle  metabolites.  We  performed  simulations  for 
the  mutant  parasite  in  the  same  way,  albeit  after  scaling  the  model 
parameters  that  govern  TCA  cycle  enzyme  activities  by  the  estimated 
increase  or  decrease  in  reaction  fluxes  of  those  enzymes.  We  used  a 
freely  available  software  package  called  BISEN  (Vanlier  et  al.,  2009)  to 
implement  the  TCA  cycle  computer  model  (Wu  et  al.,  2007). 

2.5.  Prediction  of  time-dependent  P.  falciparum  DNA  synthesis 

P.  falciparum  takes  up  nutrients  to  synthesize  DNA,  RNA,  pro¬ 
teins,  and  phospholipids.  These  processes  are  not  constitutively 


2.7.  Identifying  genes  underlying  chloroquine-induced  DNA 
synthesis  inhibition 


Chloroquine  administration  inhibits  DNA  replication  in  malaria 
parasites  (Cohen  and  Yielding,  1965;  Polet  and  Barr,  1968a).  Like¬ 
wise,  existing  transcriptome  data  (Hu  et  al.,  2010)  also  show  that 
metabolic  genes  critical  for  P.  falciparum  DNA  synthesis  are  sup¬ 
pressed  by  chloroquine  treatment.  We  sought  to  identify  the  genes 
that  are  substantially  altered  in  response  to  chloroquine  and  cause 
inhibition  of  P.  falciparum  DNA  synthesis.  Obvious  candidates  are 
genes  that  are  down-regulated  and  critical  for  P.  falciparum  DNA 
synthesis.  However,  not  all  such  genes  will  have  a  pronounced  ef¬ 
fect  on  DNA  synthesis,  which  depends  on  the  extent  to  which  an 
individual  gene  is  down-regulated  and/or  the  individual  contribu¬ 
tion  of  that  gene  to  the  DNA  synthesis  of  the  parasite.  This  suggests 
that  many  combinations  of  such  genes  may  underlie  the  metabolic 
response  associated  with  chloroquine  treatment.  Following 
Occam's  razor  (Blumer  et  al.,  1990),  we  aimed  at  identifying  the 
minimal  number  of  genes  (hereinafter  referred  to  as  the  core  genes) 
that  are  sufficient  to  completely  account  for  the  effect  of  chloro¬ 
quine  on  DNA  synthesis  by  the  parasite. 

We  identified  the  core  genes  based  on  the  Bayesian  information 
criteria  (BIC)  score  (Schwarz,  1978),  which  is  given  by, 


(„ctl— „cqn)2 


and  n  de- 


n-lnlff)  +  k'ln(n),  where  SSE  =  EiLi  maxK)_min(vnl) 
notes  the  number  of  time  points,  k  the  number  of  genes  being 


tested,  vctl  the  estimated  rate  of  DNA  synthesis  under  control 
conditions,  and  vcqn  the  estimated  rate  of  DNA  synthesis  after 
chloroquine  treatment  (with  a  =  1  for  the  k  genes  being  tested).  We 
implemented  an  algorithm  that  starts  with  a  set  of  size  k  =  1  and 
increases  the  size  of  the  test  set  as  long  as  a  given  combination  of 
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genes  for  the  new  set  (size  k  + 1)  has  a  BIC  score  lower  than  the 
minimal  BIC  score  of  the  previous  size.  We  exhaustively  repeated 
this  process  at  each  chloroquine  dose  (41,  72,  and  144  nM  for  Pf3D7, 
and  43  nM  for  PfDd2). 

3.  Results 

3.1.  Validating  the  phenotype  of  reaction  flux  estimates 

We  modified  our  previous  approach  (Fang  et  al,  2014)  to 
simulate  the  effect  of  external  stress,  as  depicted  in  Fig.  1  and 
described  in  the  Materials  and  methods  Section.  We  validated  the 
methodology  by  demonstrating  that  reaction  fluxes  obtained  in 
response  to  stress  were  sufficient  to  predict  metabolic  phenotypes 
corresponding  to  that  stress.  Fig.  2A  shows  the  biochemical  path¬ 
ways  incorporated  into  the  computer  model  used  to  predict  the 
functional  consequences  of  altered  fluxes  of  TCA  cycle  enzymes  in 
response  to  previously  studied  genetic  perturbations  (Ke  et  al., 
2015).  Fig.  2B  shows  the  predicted  alterations  in  TCA  cycle  metab¬ 
olites  in  response  to  genetic  perturbations  of  TCA  cycle  enzymes 
(Fig.  2 A,  red  circles).  We  calculated  the  altered  metabolite  con¬ 
centrations  by  scaling  the  control  parameters  in  the  TCA  cycle 
model  (Wu  et  al.,  2007)  by  the  median  of  the  predicted  changes  in 
TCA  cycle  enzyme  activities  (Fig.  2C).  Our  predictions  of  TCA  cycle 
metabolites  (Fig.  2B,  filled  bars)  qualitatively  captured  the  trends  of 
relative  change  in  experimentally  measured  TCA  cycle  metabolites 
(Fig.  2B,  open  bars).  However,  our  prediction  of  the  citrate  level 
(Fig.  2B,  filled  bars)  was  not  within  the  margin  of  error  of  the 
experimental  measurements  (Fig.  2B,  open  bars).  We  examined  all 
possible  enzymes/transporters  (e.g.,  citrate  synthase,  pyruvate  de¬ 
hydrogenase,  aspartate  aminotransferase,  citrate-malate 


antiporter)  whose  perturbation  could  yield  mitochondrial  citrate 
levels  higher  than  those  in  the  wild-type  system.  In  each  scenario, 
either  the  predicted  citrate  levels  did  not  increase  or  the  predicted 
oxaloacetate  levels  increased.  Therefore,  based  on  our  model 
analysis,  the  source  of  this  increase  in  the  experimentally  reported 
citrate  level  appears  to  be  non-mitochondrial  in  nature. 

3.2.  Chloroquine  inhibits  P.  falciparum  DNA  synthesis 

In  agreement  with  previous  experimental  observations  (Cohen 
and  Yielding,  1965;  Polet  and  Barr,  1968a),  our  simulations  sug¬ 
gested  that  chloroquine  treatment  results  in  dose-dependent  in¬ 
hibition  of  DNA  replication  (Fig.  3 A).  At  chloroquine  concentrations 
of  41  nM  (Pf3D7,  panel  A,  blue  diamonds)  and  43  nM  (PfDd2,  panel 
B,  blue  diamonds),  the  rate  of  DNA  synthesis  during  the  ring  stage 
(the  first  18  h  after  infection)  was  less  inhibited  than  it  was  during 
the  late  trophozoite  and  early  schizont  phases  (beyond  30  h  after 
infection).  These  results  are  in  agreement  with  previous  observa¬ 
tions  (Yayon  et  al.,  1983). 

Because  many  studies  suggest  that  PfDd2  is  less  sensitive  to 
chloroquine  than  Pf3D7  (Moneriz  et  al.,  2011 ),  DNA  synthesis  might 
have  been  expected  to  be  more  inhibited  in  Pf3D7  than  in  PfDd2.  In 
our  simulations,  however,  a  similar  concentration  of  chloroquine 
resulted  in  greater  DNA  synthesis  inhibition  in  PfDd2  (Fig.  3B;  blue 
diamonds  versus  black  circles,  paired-sample  t- test;  p  <  0.0001) 
than  in  Pf3D7  (Fig.  3A;  blue  diamonds  versus  black  circles,  paired- 
sample  t- test;  p  =  0.15).  These  differences  likely  stem  from  differ¬ 
ences  in  the  relative  gene  expression  ratios  of  PfDd2  and  Pf3D7. 
Notably,  the  transcriptional  profile  for  PfDd2  was  more  down- 
regulated  than  that  for  Pf3D7  at  a  similar  concentration  of  chloro¬ 
quine.  In  other  words,  the  percentage  of  down-regulated  genes 
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Fig.  1.  Schematic  representation  of  the  method  developed  to  simulate  the  effect  of  chloroquine  treatment  on  metabolic  flux  profiles  of  P.  falciparum  during  the  intraerythrocytic 
development  cycle.  The  method  shown  above  with  ‘a  =  V  is  identical  to  a  previously  proposed  method  (Fang  et  al.,  2014)  and  is  used  to  estimate  reaction  fluxes  under  control 
conditions  (black  curves  in  the  bottom  panel).  We  modified  the  previous  approach  (Fang  et  al.,  2014)  to  integrate  the  effect  of  stress  (such  as  chloroquine  treatment).  In  the 
metabolic  network  shown  above,  iq,  v2,  and  v3  are  intracellular  fluxes;  Vj  is  the  input;  ve  is  the  output;  and  vg  is  the  flux  representing  the  growth  rate.  We  simulated  the  effect  of 
external  stress  by  scaling  the  baseline  flux  of  a  reaction  (vref),  using  the  relative  change  in  the  transcriptome  data  (a)  of  that  reaction  in  response  to  stress.  The  objective  was  to  find  a 
new  intracellular  flux  distribution  (red  curves  in  the  bottom  panel)  based  on  the  modified  vref.  These  new  intracellular  fluxes  were  sufficient  to  capture  the  phenotypic  alterations 
that  occur  in  response  to  external  stress.  In  the  dotted  box,  w  represents  the  weight  of  individual  reactions,  5  represents  the  stoichiometry  of  reactions,  and  lb  and  ub  represent 
lower  and  upper  bounds  of  reactions,  respectively.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  2.  Reaction  fluxes  based  on  changes  in  gene  transcription  obtained  in  response  to  a  certain  stress  condition  are  sufficient  to  mimic  the  metabolic  phenotype  corresponding  to 
that  condition.  A)  Schematic  representation  of  the  thermodynamically  balanced  computer  model  of  mitochondrial  oxidative  phosphorylation  used  to  predict  relative  changes  in 
metabolites  on  the  basis  of  relative  changes  in  estimated  reaction  fluxes  under  wild-type  (WT)  and  knockout  (KO)  conditions;  figure  modified  and  reproduced  from  (Wu  et  al., 
2007).  The  two  red  circles  represent  the  two  genes  that  were  deleted  in  a  previous  study  (Ke  et  al.,  2015).  B)  Alterations  in  the  concentrations  of  TCA  cycle  metabolites  after 
the  deletion  of  two  TCA  cycle  enzymes.  Open  bars  show  experimentally  observed  values  and  closed  bars  the  computationally  predicted  values.  Metabolite  concentrations  were 
computed  by  using  the  model  shown  in  A  and  reaction  flux  ratios  shown  in  C.  C)  Median  relative  changes  in  TCA  cycle  enzyme  fluxes,  estimated  by  using  the  approach  shown  in 
Fig.  1,  during  the  intraerythrocytic  development  cycle  under  two  conditions  (WT  and  KO).  The  effect  of  the  stress  condition,  i.e.,  knockout  of  two  TCA  cycle  enzymes,  is  captured  by 
the  transcriptome  data  obtained  under  WT  and  KO  conditions  (Ke  et  al.,  2015).  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web 
version  of  this  article.) 
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Fig.  3.  Simulated  effect  of  chloroquine  on  rates  of  DNA  synthesis  by  the  parasite  during  the  intraerythrocytic  development  cycle  (IDC).  In  the  legend,  CQN  denotes  chloroquine.  A) 
Dose-dependent  decrease  in  the  rate  of  DNA  synthesis  for  Pf3D7.  The  DNA  synthesis  rate  at  a  CQN  dose  of  0  nM  was  estimated  from  transcriptome  data  of  Pf3D7  (Llinas  et  al.,  2006). 
The  effect  of  CQN  dose  on  the  DNA  synthesis  rate  was  simulated  by  the  approach  described  in  the  Materials  and  methods  Section,  using  transcriptome  data  of  Pf3D7  (Hu  et  al.,  2010). 
B)  The  DNA  synthesis  rate  at  a  CQN  dose  of  0  nM  was  estimated  from  transcriptome  data  of  PfDd2  (Llinas  et  al.,  2006).  The  effect  of  CQN  dose  on  the  DNA  synthesis  rate  was 
estimated  by  the  same  approach  as  that  in  A,  using  transcriptome  data  of  PfDd2  (Hu  et  al.,  2010).  For  both  strains,  the  decrease  in  DNA  synthesis  occurring  later  in  the  IDC  (beyond 
30  h  after  infection)  was  more  pronounced  than  that  occurring  early  in  the  IDC  (the  first  18  h  after  infection)  in  response  to  CQN  treatment.  To  account  for  variability  of  gene 
expression  in  our  simulations,  we  added  10%  Gaussian  noise  to  the  expression  level  of  each  gene.  The  results  were  averaged  across  twenty  independent  simulations.  Error  bars 
represent  the  standard  error  associated  with  each  prediction. 


critical  for  P.  falciparum  DNA  synthesis  was  greater  in  PfDd2  than  in 
Pf3D7  (Table  SI).  Together  with  the  highly  similar  IC50  values  re¬ 
ported  previously  (Hu  et  al.,  2010)  for  Pf3D7  (41  nM)  and  PfDd2 
(43  nM),  these  results  suggest  that  the  PfDd2  strain  used  here  may 
have  diverged  from  its  true  (chloroquine-resistant)  phenotype. 


3.3.  Genes  critical  for  P.  falciparum  DNA  synthesis 

We  found  that  81  metabolic  reactions  in  our  network  model  are 
critical  for  P.  falciparum  DNA  synthesis.  We  denoted  these  reactions 
as  critical  because  their  elimination  resulted  in  abolished  DNA 
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synthesis  (see  Section  2.6).  Of  these  81  reactions,  64  had  a  known 
gene  association  whereas  17  did  not.  These  64  metabolic  reactions 
were  encoded  by  a  total  of  44  unique  genes.  These  genes  belonged 
to  pathways  involved  in  glycolysis,  pentose  phosphate  synthesis, 
pyruvate  metabolism,  purine  metabolism,  pyrimidine  metabolism, 
redox  metabolism,  and  mitochondrial  metabolism.  In  the  online 
supplement  (see  spreadsheet),  we  provide  detailed  information, 
including  the  reaction  name,  chemical  reaction,  gene  association, 
and  pathway,  of  all  reactions  critical  for  P.  falciparum  DNA  synthesis. 

3.4.  Genes  significantly  affected  by  chloroquine  treatment 

Experimental  observations  (Cohen  and  Yielding,  1965;  Polet  and 
Barr,  1968a)  and  the  present  simulations  (Fig.  3)  indicate  that 
chloroquine  treatment  inhibits  replication  of  P.  falciparum  DNA, 
which  in  turn  suggests  that  pathways  essential  for  DNA  synthesis 
are  inhibited.  Fig.  4  shows  the  relative  change  in  the  amount  of  DNA 
synthesized  by  the  parasite  during  the  IDC  period  in  response  to 
chloroquine  treatment.  Table  S2  lists  the  actual  amounts  of  syn¬ 
thesized  DNA  for  both  strains  under  each  simulated  condition.  We 
observed  a  dose-dependent  decrease  in  the  amount  of  DNA  syn¬ 
thesized  by  the  parasite  during  the  IDC  (Fig.  4,  red  markers).  We 
used  BIC  scores  (Schwarz,  1978)  and  employed  the  algorithm 
described  in  Section  2.7  to  identify  core  genes  that  caused  this  net 
decrease  of  synthesized  DNA.  Preventing  only  the  core  genes  from 
being  down-regulated  by  chloroquine  was  sufficient  to  restore  the 
amount  of  DNA  synthesized  to  control  values  (Fig.  4,  green 
markers);  for  these  simulations,  the  reaction  expressions  of  the 
core  genes  were  set  to  be  equal  for  control  and  chloroquine-treated 
conditions.  A  list  of  the  core  genes  identified  for  each  strain  and 
chloroquine  dose  (Table  S3)  suggests  that  thioredoxin  reductase 
and  ribonucleotide  reductase  substantially  affect  the  amount  of 
synthesized  DNA.  For  PfDd2,  only  two  genes  (thioredoxin  reductase 
and  ribonucleotide  reductase)  caused  the  nearly  50%  reduction  in 
the  net  amount  of  synthesized  DNA  (Table  S3  and  Fig.  4).  These 
genes  were  not  mentioned  in  the  original  transcriptome  analysis  of 
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Fig.  4.  Simulated  effect  of  chloroquine  on  net  synthesized  DNA,  obtained  by  inte¬ 
grating  the  rate  of  DNA  synthesis  (Fig.  3)  over  the  intraerythrocytic  development  cycle. 
The  ordinate  shows  values  of  net  synthesized  DNA  (•)  relative  to  control  values  as  a 
function  of  chloroquine  dose.  The  net  amount  of  synthesized  DNA  decreases  as  the 
dose  increases  (red  symbols).  The  amount  of  synthesized  DNA  returns  back  to  control 
levels  (green  symbols)  when  a  is  set  to  1  for  the  core  genes  in  Table  S3,  while  simu¬ 
lating  the  effect  of  chloroquine  on  metabolic  fluxes.  Control  values  of  net  synthesized 
DNA  (•)  were  obtained  from  transcriptome  data  of  Pf3D7  or  PfDd2  (Llinas  et  al.,  2006) 
(Table  S2).  In  the  legend,  CTL  denotes  DNA  synthesized  under  control  conditions  after 
chloroquine  treatment,  whereas  -CG  denotes  that  after  chloroquine  treatment  with  a 
set  to  1  for  core  genes.  The  results  were  averaged  across  twenty  independent  simu¬ 
lations.  Error  bars  represent  the  standard  error  associated  with  each  prediction.  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  this  article.) 


Hu  and  colleagues,  who  used  a  fundamentally  different  approach  to 
interpret  their  transcriptome  data  (Hu  et  al.,  2010).  These  authors 
focused  on  genes  that  were  either  co-regulated  or  showed  strongly 
differential  expression  (>8  log2-fold  change)  in  response  to  a  drug. 
In  contrast,  our  methodology  of  coupling  transcriptomic  changes 
with  a  metabolic  network  model  can  predict  physiological  out¬ 
comes  as  they  relate  to  parasite  metabolism,  regardless  of  the 
magnitude  of  transcriptional  change  for  an  individual  gene.  Tran¬ 
scriptional  changes  are  not  proportional  to  changes  in  metabolic 
reaction  fluxes  (Moxley  et  al.,  2009).  Our  analysis  suggests  that  low- 
amplitude  alterations  in  these  identified  genes  (Table  S3)  were 
sufficient  to  cause  a  significant  reduction  in  parasite  DNA  synthesis. 

4.  Discussion 

4.1.  Simulating  effect  of  exogenous  stress 

Several  existing  methods  integrate  whole-genome  metabolic 
networks  and  transcriptome  data  (Colijn  et  al.,  2009;  Song  et  al., 
2014;  Wallqvist  et  al.,  2016;  Zur  et  al.,  2010).  Recently,  we  devel¬ 
oped  a  method  that  predicts  metabolic  flux  profiles  of  P.  falciparum 
during  the  IDC  (Fang  et  al.,  2014).  The  transcriptome  data  of 
P.  falciparum  obtained  under  different  stress  conditions  (e.g.,  drug 
exposure,  genetic  mutation,  etc.)  contain  information  about  its 
functional  state  and  can  be  used  to  predict  qualitative  alterations,  if 
any,  in  metabolic  flux  profiles  during  the  IDC.  However,  such  data 
are  often  collected  at  a  limited  number  of  time  points;  for  instance, 
Hu  and  colleagues  (Hu  et  al.,  2010)  obtained  P.  falciparum  tran¬ 
scriptome  data  at  only  8  selected  time  points  during  the  IDC. 
Currently,  there  is  no  approach  to  integrate  limited  data,  such  as 
those  reported  previously  (Hu  et  al.,  2010),  with  a  whole-genome 
network  to  predict  metabolic  flux  profiles  during  the  IDC.  This  is 
mainly  because  such  an  approach  requires  gene  expression  data 
throughout  the  IDC  to  predict  the  corresponding  metabolic  flux 
profiles.  Therefore,  we  developed  an  approach  that  uses  gene 
expression  data  obtained  under  a  stress  condition  to  predict  al¬ 
terations  of  metabolic  fluxes  that  correspond  to  that  condition.  We 
validated  the  developed  method  by  integrating  transcriptome  data 
obtained  under  a  given  stress  condition  and  predicting  metab- 
olomic  data  corresponding  to  that  stress  (Ke  et  al.,  2015).  Using  the 
developed  method,  we  simulated  the  effect  of  chloroquine  on  the 
metabolic  flux  profiles  of  P.  falciparum  during  the  IDC. 

4.2.  Chloroquine:  mechanism  of  action  and  resistance 

Chloroquine  inhibits  the  formation  of  hemozoin  and  leads  to 
heme  accumulation,  which  is  toxic  to  the  parasite  (Fitch,  1998). 
However,  the  specific  mechanism  underlying  this  toxicity  is  not 
known  (Ginsburg  et  al.,  1998;  Lehane  et  al.,  2012).  Some  experi¬ 
ments  suggest  inhibition  of  DNA  replication  as  a  potential  mecha¬ 
nism  of  chloroquine  action  (Cohen  and  Yielding,  1965;  Kwakye- 
Berko  and  Meshnick,  1989;  Polet  and  Barr,  1968a).  We  developed 
a  method  that  successfully  predicted  a  reduction  in  the  ability  of 
the  parasite  to  synthesize  DNA  in  response  to  chloroquine  treat¬ 
ment  (see  Figs.  3  and  4).  We  identified  twelve  distinct  genes  that 
were  substantially  affected  in  Pf3D7  and  PfDd2  after  chloroquine 
treatment;  importantly,  preventing  their  down-regulation  in 
response  to  chloroquine  treatment  was  sufficient  for  restoring  the 
P.  falciparum  DNA  synthesis  inhibited  by  chloroquine  treatment 
(see  Table  1).  These  genes  belonged  to  pathways  involved  in 
glycolysis,  redox  metabolism,  and  pyrimidine  metabolism. 

A  closer  look  at  these  genes  reveals  a  potential  mechanism  of 
chloroquine  action  that  unifies  previous  theories  and  observations 
involving  heme  (Fitch,  1998),  DNA  replication  (Polet  and  Barr, 
1968a),  and  oxidative  stress  (Lehane  et  al.,  2012;  Radfar  et  al., 
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Table  1 

Core  genes  significantly  affected  by  chloroquine. 


EC  Number3 

Gene  ID 

Enzyme  name 

Gene  abbreviation 

Pathway 

1.2.1.12 

PF14_0598 

Glyceraldehyde  3P  dehydrogenase 

GAPDH 

Glycolysis 

2. 7.2.3 

PFI1105W 

Phosphoglycerate  kinase 

PGK 

Glycolysis 

5.4.2. 1 

PF11_0208 

Phosphoglycerate  mutase 

PGM 

Glycolysis 

4.2.1.11 

PF10_0155 

Enolase 

Enolase 

Glycolysis 

1.11.1.15 

PF14_0368 

Thioredoxin  peroxidase  1 

Tpx 

Redox 

1.11.1.15 

PFL0725W 

Thioredoxin  peroxidase  2 

Tpx 

Redox 

1.17.4.1 

PF10_0154 

Ribonucleotide  reductase 

RNR 

Redox 

1.8. 1.9 

PFI1170C 

Thioredoxin  reductase 

TrxR 

Redox 

1.10.2.2 

PF14_0248 

Ubiquinol-cytochrome  c  reductase 

UQCR 

Electron  transport 

1.3.3. 1 

PFF0160C 

Dihydroorotate  dehydrogenase 

DHODH 

Electron  transport 

4.2.1. 2 

PFI1340W 

Fumarate  hydratase 

FUM 

TCA  cycleb 

2. 7.4.6 

PF13_0349 

Nucleoside  diphosphokinase 

NDPK 

Pyrimidine 

a  Enzyme  Commission  number. 
b  Tricarboxylic  acid  cycle. 


2008).  Our  proposed  mechanism,  which  is  summarized  in  Fig.  5, 
suggests  that  administration  of  chloroquine  has  two  major  meta¬ 
bolic  effects: 

i.  It  interferes  with  the  function  of  glyceraldehyde  3-phosphate 
dehydrogenase  (GAPDH),  which  reduces  the  synthesis  of  phos- 
phoenolpyruvate,  whose  carbon  fixation  activity  is  essential  for 
anaplerotic  TCA  metabolism  and  redox  balance  (Storm  et  al., 
2014). 

ii.  Its  interference  of  thioredoxin  reductase  (TrxR)  reduces  con¬ 
version  of  oxidized  thioredoxin  Trx  (S2)  to  reduced  thioredoxin 
Trx  (SH2),  which  is  essential  for  donating  electrons  to  reactions 


catalyzed  by  thioredoxin  peroxidase  (Tpx)  and  ribonucleotide 
reductase  (RNR).  Tpx  is  essential  for  converting  H2O2  (reactive) 
to  H2O  (harmless)  (Pannala  and  Dash,  2015),  and  RNR  for  con¬ 
verting  ribonucleoside  5'-diphosphates  species  (ADP,  UDP,  GDP 
and  CDP),  to  deoxy-species  (dADP,  dUDP,  dGDP  and  dCDP) 
which  are  used  for  DNA  synthesis  (Munro  and  Silva,  2012). 

The  proposed  mechanism  reveals  candidate  pathways  that  alter 
the  efficacy  of  chloroquine.  For  example,  our  analysis  points  to  TrxR 
and  GAPDH  as  candidate  molecules  that  determine  the  efficacy  of 
chloroquine.  Here,  we  focused  on  candidates  susceptible  to  a  low 
dose  of  chloroquine  because  several  pathways  will  be  inhibited  at 
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Fig.  5.  Predicted  mechanism  by  which  chloroquine  kills  the  malaria  parasite.  Chloroquine  (CQN)  leads  to  the  accumulation  of  heme.  Although  the  downstream  consequences  of  this 
effect  are  unknown,  our  analysis  suggests  that  the  accumulated  heme  leads  to  inhibition  of  glyceraldehyde  3-phosphate  dehydrogenase  (GAPDH)  and  thioredoxin  reductase  (TrxR), 
and  that  the  downstream  effects  of  these  inhibitory  events  are  the  major  determinants  of  chloroquine  efficacy.  This  hypothesis  also  provides  a  mechanistic  explanation  of  how 
chloroquine  inhibits  DNA  synthesis  (Cohen  and  Yielding,  1965;  Polet  and  Barr,  1968a).  Abbreviations:  CDP,  cytidine  diphosphate;  CO2,  carbon  dioxide;  dCDP,  deoxycytidine 
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pentose  phosphate  pathway;  RNR,  ribonucleotide  reductase;  TCA,  tricarboxylic  acid;  Tpx,  thioredoxin  peroxidase;  Trx  (S2),  oxidized  thioredoxin;  Trx  (SH2),  reduced  thioredoxin; 
UDP,  uridine  diphosphate. 
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higher  concentrations  (see  Table  S3;  chloroquine  =  144  nM).  In 
addition,  proteins  that  transport  chloroquine  out  of  food  vacuoles 
should  decrease  the  effective  concentration  of  chloroquine.  Tar¬ 
geting  these  pathways  while  administering  chloroquine  treatment 
should  enhance  the  efficacy  of  chloroquine  at  a  given  effective 
concentration. 

4.3.  Chloroquine  as  DNA  intercalator? 

The  hypothesis  that  chloroquine  kills  the  malarial  parasite  by 
DNA  intercalation  is  decades  old  (Meshnick,  1990).  Although  it  was 
the  prime  candidate  mechanism  of  chloroquine  from  the  1950s  to 
the  1970s,  it  has  now  largely  been  abandoned  for  the  past  two 
decades.  Our  model  analyses  do  not  identify  a  drug  target  or  tar¬ 
gets;  however,  they  suggest  that  the  effect  of  chloroquine  admin¬ 
istration  leads  to  inhibition  of  DNA  synthesis.  According  to  our 
model,  the  parasite  responds  to  chloroquine  by  reducing  the 
expression  of  genes  critical  for  DNA  synthesis.  Although  chloro¬ 
quine  could  specifically  intercalate  with  the  identified  core  gene  set 
in  accord  with  the  DNA  intercalation  theory  of  chloroquine 
(Meshnick,  1990;  O'Brien  et  al.,  1966a),  this  is  a  highly  unlikely 
scenario.  Heme  is  a  more  plausible  molecular  inhibitor,  given  that  it 
accumulates  in  large  amounts  within  the  parasite's  food  vacuole 
and  is  capable  of  causing  DNA  damage  (Ishikawa  et  al.,  2010; 
Rahman  et  al.,  1997),  protein  oxidation  (Kitatsuji  et  al.,  2016), 
NADPH  oxidation  (Bodaness,  1983),  and  apoptosis  (Chiabrando 
et  al.,  2014).  However,  there  is  no  direct  evidence  suggesting  that 
heme  causes  any  of  these  damages  in  P.  falciparum. 

4.4.  Implications  of  model  analyses 

Chloroquine,  which  was  once  a  highly  effective  drug,  is  now 
abandoned  in  many  malaria  endemic  regions  because  of  wide¬ 
spread  resistance.  We  identified  a  potential  mechanism  by  which 
chloroquine  inhibits  parasite  DNA  synthesis.  Knowledge  of  this 
mechanism  allowed  us  to  identify  metabolic  pathways  that  could 
be  exploited  to  increase  chloroquine  efficacy.  For  example,  our 
model  analyses  suggest  that  administration  of  H2O2  would  poten¬ 
tiate  parasite  killing  by  chloroquine,  consistent  with  in  vitro  ex¬ 
periments  (Malhotra  et  al.,  1990).  Similarly,  our  analyses  suggest 
the  hypothesis  that  targeting  parasite-specific  pathways  that  syn¬ 
thesize  phosphoenolpyruvate  and  reduce  thioredoxin  would 
potentiate  parasite  killing.  Targeting  these  pathways,  while 
administering  amodiaquine  (a  chloroquine  analogue),  should  result 
in  higher  efficacy  than  administering  chloroquine,  especially 
because  amodiaquine  accumulates  more  than  chloroquine  in  the 
parasite’s  food  vacuole  and  has  a  mechanism  similar  to  that  of 
chloroquine  (Hawley  et  al.,  1996). 

5.  Limitations  of  the  current  study 

Our  aim  was  to  identify  and  investigate  enzymes  that  contribute 
to  decreased  DNA  replication  in  response  to  chloroquine  treatment. 
Using  a  transcriptomics-based  approach,  we  can  identify  cellular 
responses  that  are  specific  to  chloroquine  by  expression  data  taken 
from  chloroquine-treated  parasites.  Although  other  antimalarial 
drugs  may  activate  similar  stress  pathways,  the  details  of  the 
response  will  be  specific  to  the  particular  conditions.  The  approach 
makes  no  assumptions  on  the  underlying  molecular  actions  that 
cause  or  affect  the  response.  The  identified  genes  are  not  likely  to  be 
affected  by  chloroquine  directly.  The  effect  of  chloroquine  on  DNA 
replication  is  likely  to  be  indirect,  being  mediated  by  heme  accu¬ 
mulation,  which  evokes  multiple  toxic  effects  in  the  cell. 

We  simulated  the  effect  of  chloroquine  over  the  entire  IDC  based 
on  the  available  transcriptome  data  taken  from  an  eight-hour 


window  during  the  trophozoite  stage.  Ideally,  for  any  specific  in¬ 
hibitor,  it  is  desirable  to  have  hourly  transcriptional  readouts  to 
predict  the  effect  of  a  systemic  perturbation  throughout  the  entire 
IDC.  In  the  case  presented  here,  the  largest  effect  of  chloroquine 
occurs  during  the  trophozoite  stage  and  the  protocol  for  collecting 
experimental  transcriptomic  data  was  designed  to  characterize 
transcriptional  changes  during  this  stage.  Our  use  of  these  data  to 
identify  changes  in  metabolite  fluxes  throughout  the  IDC  relative  to 
untreated  controls  rests  on  the  assumption  that  altering  the 
expression  of  genes  not  associated  with  any  biological  function 
specific  to  a  particular  stage  will  have  minimal  physiological  effects 
at  other  stages.  This  is  an  assumption  that  is  specific  to  P.  falciparum, 
because  malaria  relies  on  the  execution  of  a  highly  regimented 
transcriptional  program  during  the  IDC,  which  initiates  specific 
stage-dependent  cellular  functions.  These  assumptions  were  ulti¬ 
mately  corroborated  because  the  observed  metabolite  phenotype 
changes  are  mainly  associated  with  the  trophozoite  stage  and 
metabolic  alterations  during  previous  time  points  were  non¬ 
existent  or  immaterial  to  DNA  synthesis. 


6.  Conclusion 

Herein,  we  investigated  the  potential  mechanism  by  which 
chloroquine  inhibits  DNA  synthesis  in  P.  falciparum  and  kills  the 
parasite.  We  first  developed  an  approach  (that  requires  limited 
transcriptome  data)  to  simulate  the  effect  of  exogenous  stress  on 
P.  falciparum  metabolic  fluxes  during  the  IDC.  We  then  developed 
an  algorithm  employing  Bayesian  information  theory  to  identify 
potential  pathways  that  inhibit  DNA  synthesis  in  response  to 
chloroquine  treatment.  Our  analysis  suggests  that  I )  chloroquine 
inhibits  DNA  synthesis  via  inhibition  of  redox  metabolism,  carbon 
fixation,  and  anaplerotic  TCA  cycle  metabolism;  2)  P.  falciparum 
may  maintain  DNA  synthesis  after  chloroquine  treatment  by  pre¬ 
serving  a  reduced  pool  of  thioredoxin  and  by  invoking  alternate 
mechanisms  for  synthesizing  phosphoenolpyruvate,  and  3)  a  drug 
or  combination  of  drugs  that  target  redox  metabolism,  carbon  fix¬ 
ation,  and  heme  accumulation  could  be  an  effective  way  to  target 
the  parasite. 

The  method  developed  here  used  limited  transcriptome  data  to 
identify  pathways  that  underlie  the  inhibition  of  DNA  synthesis  of 
P.  falciparum  in  response  to  chloroquine  treatment.  This  suggests 
that  the  developed  method  could  also  be  suitable  to  simulate  and 
capture  the  effect  of  other  drugs  (or  stress  conditions)  on  complex 
metabolic  and  physiological  responses  of  P.  falciparum,  based  solely 
on  transcriptome  data. 
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